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Using an expansion of the order parameter over the eigen- 
functions of the hnearized first Ginzburg-Landau (GL) equa- 
tion, we obtain numerically the saddle points of the free en- 
ergy separating the stable states with different number of vor- 
tices. In contrast to known surface and geometrical barrier 
models, we find that in a wide range of magnetic fields below 
the penetration field, the saddle point state for fiux penetra- 
tion into a disk does not correspond to a vortex located nearby 
the sample boundary, but to a region of suppressed supercon- 
ductivity at the disk edge with no winding of the current, 
and which is a nucleus for the following vortex creation. The 
height of this nucleation barrier, which determines the time 
of flux penetration, is calculated for different disk radii and 
magnetic fields. 

PACS number(s): 74.24.Ha, 74.60.Ec, 73.20.Dx 



The study of magnetic flux penetration and expulsion 
in type-II superconductors has traditionally attracted 
much attention in view of important technological and 
fundamental questions concerning hysteretic behavior 
and phase transitions in bounded samples. The vortex 
creation problem is also related to phase transitions in 
superfluids. It is well known 0] that for type-II super- 
conductors (A/^ > l/-\/2, A, ^ - are the penetration and 
coherence lengths, respectively), the Meissner state be- 
comes energetically unfavorable with increasing magnetic 
field at H = Hd in comparison to the Abrikosov vortex 
lattice. In a finite system these two states, which cor- 
respond to minima of the superconductor free energy, 
are separated by a barrier. Therefore, a first-order tran- 
sition between the Meissner and Abrikosov states takes 
some time which decreases with temperature and approx- 
imately follows an Arrhenius law r cx exp{U / kT) , where 
U, T are the barrier height and sample temperature, re- 
spectively. For T = 0, the Meissner state survives up 
to the penetration field Hp and transits suddenly to the 
Abrikosov state due to dynamic instability of the order 
parameter . With decreasing the magnetic field at zero 
temperature, the vortex state remains stable up to the 
expulsion field ifg < Hd and then goes to the Meissner 
state due to vortex expulsion. 

The origin of barriers for flux penetration and expul- 
sion has been discussed during the last thirty years. Ac- 
cording to the Bean-Livingston (BL) model , the sur- 
face barrier appears due to a competition between the 
vortex attraction to the sample walls by its mirror image 
and its repulsion by screening currents. This model was 
further developed for: i) cylindrical samples, where the 



vortex shape was assumed not to be an infinite line but 
a semicircle , ii) thin disks [|| and iii) strips |^ , where 
shielding, due to finite size effects, does not decay expo- 
nentially. For samples with a non-elliptical cross section, 
the geometrical harrier arises because of Meissner screen- 
ing currents flowing on the top and bottom surfaces of 
a flat strip In addition, vortex pinning by defects 
can play an important role in the delay of vortex ex- 
pulsion or promotion of vortex penetration. It should be 
stressed, that the above mentioned barrier models, which 
are based on the London theory, do not account for the 
process of vortex formation and describe only the vortex 
motion far from the sample boundary. The Ginzburg- 
Landau (GL) theory has been previously applied for the 
study of barriers only for the ID cases of narrow wires 
and rings The approaches based on solving time- 
dependent GL equations |^ allow to treat flux penetra- 
tion (expulsion) only for magnetic fields higher (lower) 
than the penetration (expulsion) field. In this Letter, 
starting from the non-linear GL theory we present an 
approach for finding the saddle point states in thin disks 
and calculate numerically the heights of the free energy 
barriers separating the stable states with different num- 
ber of vortices. 

We consider a superconducting defect free disk with 
radius R and thickness d immersed in an insulating media 
in the presence of a perpendicular uniform magnetic field 
H. For thin disks Rd <C A^ we can neglect the distortion 
of the magnetic field, which are induced by screening and 
vortex currents, and write the GL functional as 
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where G, Gn are the free energies of the superconduct- 
ing and normal states; ^' is the complex order parameter; 
L = {~ihV — e* A/ cY /2m* is the kinetic energy operator 
for Cooper pairs of charge e* = 2e and mass m* = 2m; 
A = e^Hpjl is the vector potential of the uniform mag- 
netic field written in cylindrical coordinates 0, p; a, /3 
are the GL coefficients depending on the sample temper- 
ature. Expanding the order parameter 4" = Cii/^i 
in the orthonormal eigenfunctions of the kinetic energy 
operator L'\\}i = CiV'i p|-pd|] we go from the functional 
form (^ to the free energy written in terms of complex 
variables 

G - G„ = (a + e,)aC* + ^A%ClG]G^Gu (2) 

where the matrix elements A^^^ = J df'ip*ipj'4'kipi are cal- 
culated numerically. Note, that the sample geometry en- 
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ters in the calculations only through the eigen energies 
and eigenfunctons ipi, which are well known for the disk 
case jl^Jl^. In thin (d ^ ^) disks these eigenf unctions 
have the form 4'i={i,n} — ^xp{il4>)fn{p)^ where I is the an- 
gular momentum and the index n counts different states 
with the same I. In contrast to the approaches 
we do not restrict ourselves to the lowest Landau level 
approximation (i.e. n = 1) and expand the order param- 
eter over all eigenfunctions with e.^ < e*, where the cut- 
ting parameter e* is choosen such that increasing it does 
not influence our results. The typical number of complex 
components used are in the range N = 30 50. Thus 
the superconducting state is mapped into a 2D cluster of 
N classical particles {x,y) <-> fi?e(C), /m(C)), which is 
governed by the Hamiltonian (g) . 

To find a saddle point, which presents an extremum of 
the free energy, we use a technique similar to the eigen- 
vector following method |l3j. We start with some set of 
coefficients C . In the vicinity of this point the free energy 
5G = G{C^) — G{C) can be represented as a quadratic 
form for small deviations 5 = C" — C: 



5G — FmS„i 



DmnS*S* + C.C. 



(3) 



where F,^ = {a + e,)C™ + pA^'CjGlGi, Bmn = {a + 
im)Imn + 2/3A^"CfeC*, Dmn = f3A^/^Cf^Ci, and Imn is 
the unit matrix. The quadratic form (0), which is Hermi- 
tian, can be rewritten in normal coordinates 5m — XkQm 
as 5G = 2{-fkXk + Vkxl), where jk = Qm^m, the eigen- 
values rjk and eigenvectors are found by solving nu- 
merically the following equation 
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Moving in the direction with negative free energy gradi- 
ent ~^k we will approach a minimum of the free energy 
corresponding to the ground or a metastable state. In 
order to find a saddle point we move to a minimum of 
the free energy in all directions Xk^i = —jk/i^ + Vk) 
except one, which has the lowest eigenvalue, and for 
which we go to a maximum xi = "fi/{—e + r]i), and find 
Cm = Gm + XkQm for all k. The iteration parameter 
e > controls the convergency, which is always reached 
for any initial state close enough to a saddle point. Start- 
ing from different initial states, for which the coefficients 
C are choosen randomly, we find the saddle points for 
different magnetic fields (Fig. Due to fluctuations 
(i.e. thermal, ...) the system will be able to reach the 
saddle point and can then transfer to the other super- 
conducting state. When the magnetic field approaches 
the expulsion or penetration field, the attraction region 
of a saddle point state decreases and random searching 
becomes inefficient. Therefore, to trace the saddle point 
evolution in the vicinity of the penetration (expulsion) 
field we start from the saddle point state and increase 
(decrease) the magnetic field up to the penetration (ex- 
pulsion) field, when the lowest eigenvalue goes to zero 
(see Fig. ||, dashed curve) . 



The spatial distributions of the superconducting elec- 
tron density l^-p and velocity V = hVS - 2eA/c {'^ = 
Fexp{iS)) in the saddle point state corresponding to the 
transition from the Meissner state to the vortex state are 
depicted in Figs. ^ and ^ for different magnetic fields and 
disk radius R = 4.8^. These figures demonstrate two 
different stages in the saddle point evolution. Below the 
penetration field, the saddle point state corresponds to a 
region of suppressed superconductivity (Fig. ||(d)) with a 
minimum of l^'p, which is located at the disk boundary. 
While the minimum value of the order parameter remains 
different from zero, the vorticity L — J dldS/dl/2TT, 
where integration is performed along the disk boundary, 
equals zero and the supervelocity distributions is similar 
to that of the Meissner state (Fig. ^a)). When the order 
parameter reaches its zero value at the nucleation field 
Hn (Fig. ||, solid curve), the vorticity transits suddenly 
to L = 1. For lower magnetic fields H < Hn, the saddle 
point state presents a vortex-like state with closed veloc- 
ity circulation (Fig. ^b)). Note, that this transition is 
not followed by any discontinuity in the free energy or the 
curvature of the potential curve rj. With further decreas- 
ing the magnetic field, the saddle point corresponds to a 
vortex closer to the disk center (Fig. ^(d),(c)). This phys- 
ical picture of flux expulsion and penetration remains 
valid for other transitions L L + l with different L, in- 
dependently of the disk radius and the type (giant vortex 
or multivortex of superconducting state. 

The free energy, measured in the condensation energy 
Go = a^TrR^d/2P, is shown in Fig. ^ for the saddle point 
(dotted curves) and stable (solid curves) states for the 
disk radius R — 4.8 f . The difference between the free en- 
ergy of the saddle point state and the nearby metastable 
state corresponds to the transition barrier shown in the 
inset of Fig. js] for transition <-!■ 1. As seen from Fig. ||, 
the penetration barrier grows more slowly deep inside the 
metastable region than the expulsion barrier. Therefore, 
we expect a larger fluctuation of the penetration field 
at a finite sample temperature, which agrees with recent 
experimental observations by Geim [p^ . 

Below the nucleation field, when the saddle point state 
is similar to a vortex state, the penetration and expulsion 
barriers (see Fig. 6) can be estimated from the London 
theory, which leads to the following expression for the 
vortex free energy 

AttG. ^ TcR ' ^ R^'^A^ e 

where = a^£,'^d/2(3 = Gq^'^/ttR'^, p is the radial 
vortex position, rc ~ ^ is the vortex core radius, $ = 
ttHR^/^o is the unitless magnetic flux, and $o = hc/2e 
is the flux quantum. Note that: i) the expulsion field 
He — ^o/ttR^, ii) the vortex position ps = Ry/l — 1/$ 
in the saddle point, and iii) the BL expulsion barrier 
U = Gi(p,)-Gi(0) = 47rG^($-l-?n$) do not depend 
on the vortex core energy, which is represented by the 
first term in Eq. (0). As seen from Fig. @(a), the London 
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theory prediction for the expulsion barriers are confirmed 
by the GL theory in the hmit of large disks ^. We 
extended the BL model to arbitrary R/^ by taking into 
account the spatial nonuniformity of the modulus of the 
order parameter, which obeys the first GL equation 

with the boundary condition {dF/dp)p=ji = 0. Follow- 
ing the BL model we assume that the phase distribution 
is created by a vortex and its mirror image, which are 
located at the distances pv,R^/pv from the disk center 
1^,^, respectively. Solving Eq. ^ numerically for dif- 
ferent vortex positions py we find the expulsion (solid 
circles) and penetration (open circles) barriers shown in 
the inset of Fig. ^. Below the nucleation field there is an 
excellent quantative agreement between our GL theory 
and this improved BL model. Nevertheless, this model 
breaks down in the range 7J„ < H < Hp. Note, that 
the barrier height at H = Hn increases with disk radius 
and the role of the nucleation barrier may become even 
more important in macroscopic systems, where possible 
3D (for d> £_) and demagnetization (for Rd > A^) effects 
must also be taken into account. In the unitless variables 
2ttH^R/^q (<i>/<i>o), the penetration (expulsion) barriers 
measured in Gj, are proportional to the disk thickness 
and increase slightly with the disk radius. 

In conclusion, we have demonstrated that in a wide 
range of magnetic fields i?„ < H < Hp the saddle point 
state presents a vortex nucleus, which is a region of sup- 
pressed superconductivity surrounded by a background 
of Meissner state, which transits to a vortex state at 
H < Hn- We have found the penetration field and 
the corresponding nucleation barriers for thin disks. For 
lower magnetic fields Hg < H < Hn, the saddle point 
state can be reasonably described by the conventional 
London theory. We also extended the BL model to finite 
disk radius. 
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FIG. 1. Schematical view of the free energy in functional 
space depicting two minima with L = 2 and L = 3 vortices 
and the saddle point connecting them. The Cooper pair den- 
sity of these three states are shown in the insets. 

FIG. 2. The minimal density of superconducting electrons 
|\['|^ at a point of the disk boundary (solid curve), measured 
in —a/ (5, and the curvature 77 (dashed curve) of the free en- 
ergy curve for the transition between the Meissner and vortex 
states. 
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FIG. 3. Density of superconducting electrons l*!'!^ in the 
saddle point states corresponding to the transition 0^1 
for the disk radius R = 4.8^ and different magnetic fields. 
The signs (-I-), (-) denote the maxima and minima of 
respectively. 
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FIG. 4. The velocity of superconducting electrons in the 
saddle point states corresponding to the transition <-» 1 for 
the disk radius R = 4.8^ and for two different magnetic fields. 

FIG. 5. The free energy of the stable (solid) and sad- 
dle point (dotted curves) states as function of the magnetic 
field, measured in Hc2 ~ ^o/2'k(,^, for different angular mo- 
menta L. The inset shows the expulsion (solid circles) and 
penetration (open circles) barriers found with our modified 
Bean-Livinston model (symbols) which are compared to our 
GL approach (curves) for the transition <-» 1. 

FIG. 6. The expulsion (a) and penetration (b) barrier as 

function of the magnetic field for different disk radii. The 
dashed curve corresponds to the London limit for the expul- 
sion barrier. The solid circles correspond to the penetration 
barrier at the nucleation field. 
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